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'"IT' ' Abstract 

Y^- A complete analysis of multi-mode bosonic Gaussian channels is proposed. 

j-H , We clarify the structure of unitary dilations of general Gaussian channels in- 

C^ I volving any number of bosonic modes and present a normal form. The maxi- 

mum number of auxiliary modes that is needed is identified, including all rank 
deficient cases, and the specific role of additive classical noise is highlighted. 
By using this analysis, we derive a canonical matrix form of the noisy evolution 
^ ' of n-mode bosonic Gaussian channels and of their weak complementary coun- 

terparts, based on a recent generalization of the normal mode decomposition 
for non-symmetric or locality constrained situations. It allows us to simplify 
O ' the weak-degradability classification. Moreover, we investigate the structure 

Tjj- . of some singular multi-mode channels, like the additive classical noise chan- 

O I nel that can be used to decompose a noisy channel in terms of a less noisy 

one in order to find new sets of maps with zero quantum capacity. Finally, the 
two-mode case is analyzed in detail. By exploiting the composition rules of 
two-mode maps and the fact that anti-degradable channels cannot be used to 
rN I transfer quantum information, we identify sets of two-mode bosonic channels 

j^ ■ with zero capacity. 
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Bosonic Gaussian channels are ubiquitous in physics. They arise whenever a har- 
monic system interacts linearly with a number of bosonic modes which are inaccessi- 
ble in principle or in practice |[Il|2l[3l|4l[51[6l|71|. They provide realistic noise models 
for a variety of quantum optical and solid state systems when treated as open quan- 
tum systems, including models for wave guides and quantum condensates. They play 
a fundamental role in characterizing the efficiency of a variety of tasks in continuous- 
variables quantum information processing [ 8 1 , including quantum communication |l9l 
and cryptography [TTOll . Most importantly, communication channels such as optical 
fibers can to a good approximation be described by Gaussian quantum channels. 

Not very surprisingly in the light of the central status of such quantum channels, 
a lot of effort has been recently devoted to studying their properties (see Ref . [4J for 
a review), based on a long tradition of work on Gaussian channels (Hill [3. Specif- 
ically, from a quantum information perspective, a key question is whether or not 
a channels allows for the reliable transmission of classical or quantum information 
llallSIIIlinilllllllBlIiainilll. significant progress has been made in this re- 
spect in recent years, although for some important cases, like the thermal noise chan- 
nel modelling a realistic fiber with offset noise, the quantum capacity is still not yet 
known. In this context, the degradability properties represent a powerful tool to sim- 
plify the quantum capacity issue of such Gaussian channels. Indeed, in Refs. [[T6llT7l 
it has been shown that with some (important) exceptions, Gaussian channels which 
operate on a single bosonic mode (i.e., one-mode Gaussian channels) can be classi- 
fied as weakly degradable or anti-degradable. This paved the way for the solution of 



the quantum capacity fTO'l for a large class of these maps fT5\. 

Here, first we propose a general construction of unitary dilations of multi-mode 
quantum channels, including all rank-deficient cases. We characterize the minimal 
noise maps involving only true quantum noise. Then, by using a generalized nor- 
mal mode decomposition recently introduced in Ref. ||20l , we generalize the results 
of Refs. |fT6l [m concerning Gaussian weak complementary channels to the multi- 
mode case giving a simple weak-degradability/anti-degradability condition for such 
channels. The paper ends with a detailed analysis of the two-mode case. This is 
important since any n-mode channel can always be reduced to single-mode and two- 
mode components [|20l . We detalize the degradability analysis and investigate a use- 
ful decomposition of a channel with the additive classical noise map that allows us to 
find new sets of channels with zero quantum capacity. 

1 Multi-mode bosonic Gaussian channels 

Gaussian channels arise from linear dynamics of open bosonic system interacting 
with a Gaussian environment via quadratic Hamiltonians. Loosely speaking, they 
can be characterized as CPT maps that transform Gaussian states into Gaussian states 

1.1 Notation and preliminaries 

Consider a system composed by n bosonic modes having canonical coordinates 
Qi, Pi,- ■ ■ ,Qn,Pn- The canonical commutation relations of the canonical coordi- 
nates, [Rj, Rj'] = i{a2n)j,j', where R := {Qi,- ■ ■ ,Qn; Pi,- ■ ■ , Pn), are grasped by 
the 2n x 2n commutation matrix 



(^2n 



K 

-Un 



(1) 



when this order of canonical coordinates is chosen, (here U^ is the n x n identity 
matrix) [|3l IH |2T1- Even though different reordering of the elements of R will not 
affect the definitions that follow, we find it useful to assume a specific form for cr2„. 
One defines the group of real symplectic matrices Sp{2n, H) as the set of 2n x 2n 
real matrices S which satisfy the condition 

Sa2nS = a2n ■ (2) 

Since Det[cr2n] = 1. and a2n = —<^2n, any symplectic matrix S has Det[S'] = 1 
and it is invertible with S^^ G Sp{2n,I{). Similarly, one has 5^ G 5^(2^,11). 
Symplectic matrices play a key role in the characterization of bosonic systems. In- 
deed, define the Weyl (displacement) operators as V{z) = V^^i^—z) := exp[iRz] 
with z := (xi, X2,- ■ ■ , x„, yi, ?/2, ■ ■ ■ , 2/n)^ being a column vector of H^"-. Then it 
is possible to show [1] that for any S G Sp{2n, R) there exists a canonical unitary 
transformation U which maps the canonical observables of the system into a linear 
combination of the operators Rj, satisfying the condition 

f/t v(z) U = V{Sz) , (3) 



for all z. This is often referred to as metaplectic representation. Conversely, one can 
show that any unitary IJ which transforms V^z) as in Eq. ([3]) must correspond to an 
Se ^p(2n,Il). 

Weyl operators allow one to rewrite the canonical commutation relations as 

V{z)V{z') = expi-^z'^a^nz'jViz + z') , (4) 

and permit a complete descriptions of the system in terms of (characteristic) com- 
plex functions. Specifically, any trace-class operator (in particular, any density 
operator) can be expressed as 

(^'^(0'-)^-^), (5) 

where d'^"^z := dxi ■ ■ ■ dxndyi ■ ■ ■ di/n and 0(0; z) is the characteristic function asso- 
ciated with the operator defined by 

4>{Q;z):=Tr[QV{z)]. (6) 

Within this framework a density operator p of the n modes is said to represent a 
Gaussian state if its characteristic function (hip; z) has a Gaussian form, i.e.. 



0(/); z) = exp[-|Z 'yz + im z] , (7) 

with m being a real vector of mean values rrij := Tr:[ pRj], and the 2n x 2n real 
symmetric matrix 7 being the covariance matrix [HI SlEl of p. For generic density 
operators p (not only the Gaussian ones) the latter is defined as the variance of the 
canonical coordinates R, i.e.. 



7j,i' : = Tr p { (/?j - mj ) , (i?jv - mjv ) } 



(8) 



with {■, ■} being the anti-commutator, and it is bound to satisfy the uncertainty rela- 
tions 

7 ^ ^(^2n , (9) 

with a2n being the commutation matrix ([T]). Up to an arbitrary vector m, the uncer- 
tainty inequality presented above uniquely characterizes the set of Gaussian states, 
i.e. any 7 satisfying ^ defines a Gaussian state. Let us first notice that if 7 satis- 
fies © then it must be (strictly) positive definite 7 > 0, and have Det[7] ^ 1. From 
Williamson theorem [|23l it follows that there exists a symplectic S E Sp{2n, R) 
such that 



7 = ^ 



D 
D 



S^ , (10) 



where D := diag((ii, ■ ■ ■ , dn) is a diagonal matrix formed by the symplectic eigen- 
values dj ^ 1 of 7. For S = t2n Eq. (fTOl) gives the covariance matrix associated with 



thermal bosonic states. This also shows that any covariance matrix 7 satisfying Q 
can be written as 

7 = ^^^ + A, (11) 

with A ^ Qj. The extremal solutions of Eq. ([TT]) . i.e., 7 = SS'^, are minimal 
uncertainty solutions and correspond to the pure Gaussian states of n modes (e.g., 
multi-mode squeezed vacuum states). They are uniquely determined by the condition 
Det[7] = 1 and satisfy the condition [flSlI 

7 = -0-2n(7~^)0"2n- (13) 

1.2 Bosonic Gaussian channels 

In the Schrodinger picture evolution is described by applying the transformation 
to the states (i.e., the density operators), p ^^ ^{p)- In the Heisenberg picture 
the transformation is applied to the observables of the system, while leaving the 
states unchanged, i— ^ $//(©). The two pictures are related through the identity 
Tr[$(p)e] = Tr[p$i^(e)], which holds for all p and 0. The map $// is called the 
Jaa/of $. 

Due to the representation (|5]) and Q any completely positive, trace preserving 
(CPT) transformation on the n-modes can be characterized by its action on the Weyl 
operators of the system in the Heisenberg picture (e.g., see Ref. [Tf|). In particular, 
a bosonic Gaussian channel (BGC) is defined as a map which, for all z, operates on 
V{z) according to Q 

V{z) I — > ^H{y{z)) ■= V{Xz) exp[-lz'^Yz + iv^z] , (14) 

with V being some fixed real vector of R^", and with Y,Xe ]R2nx2n \)Q[^g some 
fixed real 2n x 2n matrices satisfying the complete positivity condition 

y ^ iS with S := a2n - X^a2nX . (15) 

In the context of BGCs the above inequality is the quantum channel counterpart of 
the uncertainty relation ^. Indeed up to a vector v, Eq. (fTSi) uniquely determines the 
set of BGCs and bounds Y to be positive-semidefinite, Y ^ 0. However, differently 
from ^ in this case strict positivity is not a necessary prerequisite for Y. A com- 
pletely positive map defined by Eqs. (fT4l) and ([T?] ) will be referred to as bosonic Gaus- 
sian channel (BGC). As mentioned before, such a map is a model for a wide class 
of physical situations, ranging from communication channels such as optical fibers, 
to open quantum systems, and to dynamics in harmonic lattice systems. Whenever 
one has only partial access to the dynamics of a system that can be well-described 



'This is indeed the matrix 



D-t, 



S^ (12) 



with D as in Eq. (flOb which is positive since D ^ in- 



by a time evolution governed by a Hamiltonian that is a quadratic polynomial in the 
canonical coordinates, one will arrive at a model described by Eqs. (fT4l) and (fTSHi. 

An important subset of the BGCs is given by set of Gaussian unitary transforma- 
tions which have Y = 0, X e Sp{2n, R), and v arbitrary. They include the canonical 
transformations of Eq. ([3]) (characterized by t> = 0), and the displacement transfor- 
mations (characterized by having X = t2n and v arbitrary). The latter simply adds 
a phase to the Weyl operators and correspond to unitary transformations of the form 
^h{V{z)) := V{-v)V{z)V{v) = V{z) expliv'^z]. 

In the Schrodinger picture the BGC transformation (fT4)) induces a mapping of the 
characteristic functions of the form 

0(p; z) ^^ 0($(/5); z) := 0(p; Xz) expl-^z'^Yz + iv"" z\ , (16) 

which in turn yields the following transformation of the mean and the covariance 
matrix 

m I — > X'^m + V , 
7 I — > X^7X + F. (17) 

Clearly, BGCs always map Gaussian input states into Gaussian output states. 

For purposes of assessing quantum or classical information capacities, output en- 
tropies, or studying degradability or anti-degradability of a channel [|T6l[T7iri4llT5l , 
the full knowledge of the channel is not required: Transforming the input or the out- 
put with any unitary operation (say, Gaussian unitaries) will not alter any of these 
quantities. It is then convenient to take advantage of this freedom to simplify the 
description of the BGCs. To do so we first notice that the set of Gaussian maps is 
closed under composition. Consider then $' and $" two BGCs described respec- 
tively by the elements X', F' , v' and X", Y" , v" . The composition $" o $' where, in 
Schrodinger representation, we first operate with $' and then with $", is still a BGC 
and it is characterized by the parameters 

V = {X"fv' + v\ 

X = x'x", 

Y = {X"fY'X" + Y". (18) 

Exploiting these composition rules it is then easy to verify that the vector v can 
always be compensated by properly displacing either the input state or the output 
state (or both) of the channel. For instance by taking X" = t2n, Y" = and 
v" = —v', Eq. (fTSi) shows that $' is unitarily equivalent to the Gaussian channel $ 
which has v = and X = X', Y = Y' . Therefore, without loss of generality, in the 
following we will focus on BGCs having f = 0. 

More generally consider the case where we cascade a generic BGC $' described 
by matrices X', Y' as in Eq. (fT5l) with a couple of canonical unitary transformation 
Ui and U2 described by the symplectic matrices 5*1 and 5*2 respectively. The resulting 



^This set does not contain ideal Gaussian measurements, like optical homodyning |[3T1 . 



BGC $ is then described by the matrices 



X = ^i(X')^2, (19) 

Y = S'i{Y')S2. 



For single mode (n = 1) this procedure induces a simplified canonical form 11131 
[TTl [T4ll which, up to a Gaussian unitarily equivalence, allows one to focus only on 
transformations characterized by X and Y which, apart from some special cases, are 
proportional to the identity. In this paper we will generalize some of these results to 
an arbitrary number of modes n. To achieve this goal, in the following section we 
first present an explicit dilation representation in which the mapping ((14)) is described 
as a (canonical) unitary coupling between the n modes of the system and some extra 
environmental modes which are initially prepared into a Gaussian state. Then we will 
introduce the notion of minimal noise channel, showing a useful decomposition rule. 



2 Unitary dilation theorem 

In this section we introduce a general construction of unitary dilations of multi-mode 
quantum channels. Specifically we show that a CPT channel acting on n modes is a 
BGC if and only if it can be realized by invoking i ^ 2n additional (environmental) 
modes E through the expression 

<l>(p)=Tr^[f/(/)®/)E)f/t], (20) 

where p is the input n-mode state of the system, pe is a Gaussian state of an envi- 
ronment, [/ is a canonical unitary transformation which couples the system with the 
environment, and Trg denotes the partial trace over E. In case in which pE is pure, 
Eq. (|20|) corresponds to a Stinespring dilation [ J24l of the channel $, otherwise it is a 
physical representation analogous to those employed in Refs. |[T6l [T71 for the single 
mode case. 



2.1 General dilations 

In this subsection, we will construct Gaussian dilations, including a discussion of all 
rank-deficient cases, and will later focus on dilations involving the minimal number 
of modes. To proceed, we will first establish some conventions and notation. To start 
with, we write the commutation matrix of our n + i modes in the block structure 



o- := a2n © o-2i 



0-2„ 

erg 



^^"^ (21) 



where a2n and a^ are 2n x 2n and 2i x 2i commutation matrices associated with 
the system and environmental modes, respectively. For a2n we assume the structure 
as defined in Eq. ([T]). For a^, in contrast, we do not make any assumption at this 



point, leaving open the possibility of defining it later or|j. Accordingly, the canonical 
unitary transformation U of Eq. (l20l) will be uniquely determined by a 2(?2 + £) x 
2(n + (i) real matrix S E Sp{2{n + i),R) of block form 



S :-- 



S3 S4 



(22) 



which satisfies the condition 



SaS 



T 



a . 



E ^T 



Si 0'2n ■§! + S2 O^^ S2 — 0'2n 



■Si 0"2n S 



T 



S2 0"2£ "^4 



E ,.T 



(23) 



■S3 '^2n ■Ss + ■S4 <^2e. ■^4 ~ '^: 



2£ ' 



In the above expressions, Si and S4 are 2n x 2n and 2£ x 2£ real square matrices, 
while S2 and s|^ are 2n x 2£ real rectangular matrices. Introducing then the covariance 
matrices 7 ^ 2(j2„ and 7^ ^ icr^ of the states p and p^;, the identity (|20l ) can be 
written as 



5 



7 

7£; 



S^ 



Si 7 sf + S2 7e -s^ 



X'^X + Y, 



(24) 



2n 



where |2n denotes the upper principle submatrix of degree 2n, and where X,Y E 
j^2nx2?i satisfying the condition (fT5l) are the matrices associated with the channel $. 
In writing Eq. (I24l) we use the fact that due to the definition (|2T]) the covariance matrix 
of the composite state p® pe can be expressed as 7 © 7^;. With these definitions, the 
first part of the unitary dilation property (|20l) can be written as follows: 

Proposition 1 (Unitary dilations of Gaussian channels) Let ^e be the covariance 
matrix of a Gaussian state of t modes and let S G Sp{2{n + £), R) be a symplectic 
transformation. Then there exists a symmetric 2n x 2n-matrix Y ^ and a2nx 2n- 
matrix X satisfying the condition d75l) . such that Eq. ^24\ holds for all 7. 

Proof: The proof is straightforward: We write S in the block form (|22)) and take 
X = si and Y = S27£;.s|'. Since 7^; is a covariance matrix of I modes, 7^; — ia2i ^ 
and therefore S2{'1e — ic^dsl ^ 0. This leads to Eq. (flSl) through the identity the 
symplectic condition Sia2ns\ + S2(J2is^ = cr2n which follows by comparing the up- 
per principle submatrices of degree n of both terms of Eq. (|23l) . ■ 



This proves that any CPT map obtained by coupling the n modes with a Gaussian 
state of £ environmental bosonic modes through a Gaussian unitary t/ is a BGC. The 
converse property is more demanding. In order to present it we find it useful to state 
first the following 



^With this choice the canonical commutation relations of the n + £ mode read as [Rj , Rj'] — ic 



where R := (Qi,--- ,Q«;^i,--- ,-Pn;^i, 
the j-th system mode and with and fi, ■ • 

have^:= (Qi,--- ,Q„;A,--- ,Pn]Qf, 



Qf, Pi\ ■ • ■ ; Qf ^ Pf of the environmental modes. For instance, taking ct^ 



f2i) with Qj, Pj being the canonical coordinates of 
r2e being some ordering of the canonical coordinates 

(72£ corresponds to 



\E. pE pE\ 



Lemma 1 (Extensions of symplectic forms) Let, for some skew symmetric a 2^, si 
and S2 be 2n x 2n and 2n x 2C, real matrices forming a symplectic system, i.e., 
si o"2n sf + S2 cr^i s^ = o"2n- Then we can always find real matrices S3 and S4 such 
that S ofEq. d22l) is symplectic with respect to the commutation matrix (|27]). 

Proof: Since the rows of S form a symplectic basis, given si and S2 (an incomplete 
symplectic basis), it is always possible to find S3 and S4 as above. The proof easily 
follows from a skew-symmetric version of the Gram-Schmidt process to construct a 
symplectic basis ll25l . For a special subset of BGCs, in Sec. 12.51 we will present an 
explicit expression for S based on a simplified (canonical) representation of the X 
matrix that defines $. See also Appendix lAl ■ 



Due to the above result, the possibility of realizing unitary dilation Eq. (1201) for a 
generic BGC described by the matrices X and F ^ iS = i(o"2„ — X'^a2nX), can be 
proven by simply taking si = X'^ and finding some 2n x 2i. real matrix S2 and an 
^-mode covariance matrix 7^; ^ za^ that solve the equations 

S2 (T21 sl = Cr2n - Si a2n sJ = E, (25) 

S2 7i?4 = y- (26) 



With this choice in fact Eq. (1241) is trivially satisfied for all 7, while si and S2 can 
be completed to a symplectic matrix S E Sp{2{n + £),]R). The unitary dilation 
property (|20l ) can hence be restated as follows: 

Theorem 1 (Unitary dilations of Gaussian channels: Converse implication) For 

any real 2n x 2n-matrices X and Y satisfying the condition ( 1731) . there exist I smaller 
than or equal to 2n, S G Sp{2{n + £), R), and a covariance matrix ■je ofi modes, 
such that Eq. i[24\) is satisfied. 

Proof: As already noticed the whole problem can be solved by assuming Si = 
X'^ and finding S2 and ^e that satisfy Eqs. (|25l ) and (|26l ). We start by observing 
that the 2n x 2n matrix S defined in Eq. (fT5l ) is skew-symmetric, i.e., S = — S^. 
Moreover according to Eq. (fT5l) its support must be contained in the support of Y, 
i.e., Supp[S] C Supp[F]. Consequently given k := rank[F] and r := rank[S] as 
the ranks of Y and S, respectively, one has that k ^ r. We can hence identify three 
different regimes: 

(i) k = 2n, r = 2n, i.e., both Y and S are full rank and hence invertible. Loosely 
speaking, this means that all the noise components in the channel are basically 
quantum (although may involve classical noise as well). 

(ii) k = 2n and r < 2n, i.e., Y is full rank and hence invertible, while S is singular. 
This means that the some of the noise components can be purely classical, but 
still nondegenerate. 

(iii) 2n > k ^ r, i.e., both Y and S are singular. There are noise components with 
zero variance. 



Even though (i) and (ii) admit similar solutions, it is instructive to analyze them 
separately. In the former case in fact there is a simple direct way of constructing a 
physical dilation of the channel with i = n environmental modes. 

(i) Since S is skew-symmetric and invertible there exists an O G 0{2n, H) or- 
thogonal such that 



oj:o^ 



/i 
-/i 



(27) 
n (see page 107 in Ref. 

(28) 
(29) 



which corresponds to Eq. (|25l) for i = n. Since si = X'^, Lemma 1 guarantees that 
this is sufficient to prove the existence of S. The condition (|24)) finally follows by 
taking ^e = KYK^ which is strictly positive (indeed K is invertible and F > 
because it has full rank) and which satisfies the uncertainty relation ^, i.e., 



where ^ = diag(/ii, ■ ■ ■ , /i„) and fit > for alH = 1, 
^M ). Hence K := M'^'^O with M := /x © // satisfies 



Taking then S2 := K ^ we ge 



S2 0^2n sl = K ^02nK ^ 



Y >%L 



lE 



KYK^ ^ iKT.K^ 



lCr2r. 



(30) 



This shows that the channel admits a unitary dilation of the form as specified in 
Eq. (|20l) with i. = n environmental modes with commutation matrix, a2n = o"2n - see 
discussion after Eq. (|2TI) . Such a solution, however, will involve a pure state pE only 

ifDet[7s] = l,i.e.. 



Det[F]Det[ir]^ 



1 



Det[r] = Det[S] . 



(31) 



When Det[7£;] > 1, i.e., Det[F] > Det[S], we can still construct a pure dilation by 
simply adding further n modes which purify the state associated with the co variance 
matrix 7^; and by extending the unitary operator (j associated with S as the identity 
operator on them. For details see the discussion of case (ii) given below. This corre- 
sponds to constructing a unitary dilation (|20l) with the pure state pE being defined on 
i = 2n modes. 



(ii) In this case Y is still invertible, while S is not. Differently from the approach 
we adopted in solving case (i), we here derive directly a Stinespring unitary dila- 
tion, i.e., we construct a solution with a pure 7^ that involves i = 2n environmental 
modes. In the next section, however, we will show that, dropping the purity require- 
ment, one can construct unitary dilation that involves pe with only i = 2n — r/2 
modes. 



"^From now on, the symbol A ^ will be used to indicate the transpose of the inverse of the matrix 

A,i.e.,A-^ -.^{A-Y^iA^y'- 



lO 



To find S2 and 7^ which solve Eqs. (1251) and (1261) . it is useful to first transform Y 
into a simpler form by a congruent transformation, i.e., 



CYC^ = t2n 



(32) 



with C e G/(2n, R) being not singular, e.g., C := F"^/^. From Eq. ([B]) it then 
follows that 

a2„^^S', (33) 

with S' := y^i/s j]y-i/2 being skew-symmetric (i.e., E' = — (S')^) and singular 
with rank[S'] = rank[S] = r Il26ll . We then observe that introducing 



„ _ vl/2 „/ 

the conditions (|25l ) and (|26l ) can be written as 

4 7s (4)^ = 



11 2n 



(34) 



(35) 
(36) 



Finding S2 and 7^; which satisfy these expressions will provide us also a solution of 
Eqs. ([25]) and (|26l). 

As in the case of Eq. (l27l) . there exists an orthogonal matrix O G 0{2n, R) 
which transforms the skew- symmetric matrix S' in a simplified block form. In this 
case however, since S' is singular, we find ff26l 



0X0 



irxT 






/i 





}r/2 








}n-r/2 


-/i 








}r/2 








}n-r/2 



(37) 



where now /i = diag(yUi, ■ ■ ■ , jJir/2) is the r/2 x r/2 diagonal matrix formed by the 
strictly positive eigenvalues of |S'| which satisfy the conditions 1 ^ /ij > 0, this 
being equivalent with 

K/2 > /^, (38) 

as a consequence of inequality (|33] ). Define then K := M~^/^ O with 



M 



^ 











fln-r/2 





/i 








ln-r/2 



It satisfies the identity 



KE'K 



iT^T 






]lr/2 











-trl2 















}r/2 
}n — r/2 
}r/2 
}n — r/2. 



}r/2 
}n — r/2 
}r/2 
} n — r/2. 



(39) 



(40) 



11 



To show that Eqs. (l35l) and (|36l) admit a solution we take i = 2n and write a^ = 
c"2n ©cr2n = '74n with a2n as in Eq. ([T])- With these definitions the 2n x 4n rectangular 
matrix S2 can be chosen to have the block structure 



with A being the following 2n x 2n symmetric matrix 



(41) 



A = A^ 















ln-r/2 














fln-r/2 



}r/2 

} n — r/2 

}r/2 

}n — r/2. 



(42) 



By direct substitution one can easily verify that Eq. (1351) is indeed satisfied, see Ap- 
pendix |B] for details. Inserting Eq. (|4T1) into Eq. (l36l) yields now the following equa- 
tion 



a + A6^ + 6 A^ + A (3 A^ = M-^ 
for the 4n x 4?t, covariance matrix 



7e 



a (5 
<5^ P 



(43) 



(44) 



see Appendix O for details. A solution can be easily derived by taking 



a = P 



\ "'' 













^ln-r/2 







/^-^ 








^i-n-r/2 



with C = 5/4 and 






/(/i-^) 








/(an-./2) 


fif^-') 











/(ean-r/2) 



}r/2 




}n-r/2 




}r/2 




}n-r/2 




1 }r/2 




}r/2 


V2 




} n - r/2 



(45) 



(46) 



with f{6) := —(6'^ — 1)^/^. For all diagonal matrices /i compatible with the con- 
straint (l38l) the resulting 7^; satisfies the uncertainty relation 'Je ^ ^cr4„. Moreover 
since it has Det[7£;] = 1, this is also a minimal uncertainty state, i.e., a pure Gaussian 
state of 2n modes. It is worth stressing that for r = 2n, i.e., when also the rank of 
S is maximum, the above solution provides an alternative derivation of the unitary 
dilation discussed in the part (i) of the theorem. In this case the covariance matrix 7^ 
has block elements 



a = P 



/i 



-1 







fi' 



}n 


6 = 


/(/i-i) " 


}n 


}n ' 


/(/^'^) 


}n 



(47) 
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where /i is now a strictly positive n x n matrix, while Eqs. (1341) and (1411) yield 



s, := yV^oT 



r /iV2 














^V2 









} n 
} n. 



(48) 



(iii) Here both Y and S are singular. This case is very similar to case (ii). Here, 
the dilation can be constructed by introducing a strictly positive matrix F > which 
satisfies the condition 



UYU = Y. 



(49) 



with n being the projector onto the support of Y. Such a Y always exists (Y = 
Y + (t — n)). By construction, it satisfies the inequality Y ^ Y ^ iE. According 
to Sec. ll.2[ Y and X define thus a BGC. Moreover, since Y is strictly positive, it has 
full rank. Therefore, we can use part (ii) of the proof to find a2n x 2i matrix S2 and 
7e ^ 'i'0-2i which satisfy the conditions (|25l) and (|26l ). i.e. 



S2 lE sl 



Y . 



A unitary dilation for the channel F, X is then obtained by choosing 7^; 
S2 = ns2. In fact from Eq. (ISTI) we get 



S2 7S 52^ 



ns2 7£;s^n 



nFn = F. 



while from Eq. (|50l) 

S2 af^ s^ = n S2 af^ s^ n = n s n = s ^ 

where we have used the fact that Supp[S] C Supp[F]. ■ 



(50) 
(51) 

^E and 



(52) 



(53) 



In proving the second part of the unitary dilations theorem we provided explicit 
expressions for the environmental state pe of Eq. (|20l) . Specifically such a state is 
given by the pure 2n mode Gaussian state pe characterized by the covariance matrix 
7e of elements (|45l) and (|46|) . A trivial observation is that this can always be replaced 



by the 2n modes vacuum state |0)(0| having the covariance matrix 7^^ = t2n- 
This is a consequence of the obvious property that according to Eq. (fTTI) all pure 
Gaussian states are equivalent to |0)(0| up to a Gaussian unitary transformation. 
On the level of covariance matrices, Gaussian unitaries correspond to symplectic 
transformations. For a remark on unitarily equivalent dilations, see also Appendix D. 
Hence, by means of a congruence with an appropriate symplectic transformation, we 
immediately arrive at the following corollary: 

Corollary 1 (Gaussian channels with pure Gaussian dilations) Any n-mode Gaus- 
sian channel $ admits a Gaussian unitary dilation d20l) with pE = |0)(0| being the 
vacuum state on 2n modes. 
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2.2 Reducing the number of environmental modes 

An interesting question is the characterization of the minimal number of environmen- 
tal modes i that need to be involved in the unitary dilation (|20l) . From Theorem 1 we 
know that such number is certainly smaller than or equal to twice the number n of 
modes on which the BGC is operating: We have in fact explicitly constructed one of 
such representations that involves i = 2n modes in a minimal uncertainty, i.e., pure 
Gaussian state. We also know, however, that there are situationgjin which i can be re- 
duced to just n: This happens for instance for BGCs $ with rank[F] = rank[S] = 2n, 
i.e., case (i) of Theorem 1 . In this case one can represent the channel $ in terms of a 
Gaussian unitary coupling with i = n environmental modes which are prepared into 
a Gaussian state with covariance matrix 

'yE = KYK^, (54) 

- see Eq. (|30l ). In general, this will not be of Stinespring form (not be a pure unitary 
dilation) since 'Je is not a minimal uncertainty covariance matrix. In fact, forn = 1 
this corresponds to the physical representation of $ of Refs. [17]. However if Y and 
X satisfy the condition (|3TI) . our analysis provides a unitary dilation involving merely 
i = n modes in a pure Gaussian state. 

We can then formulate a necessary and sufficient condition for the channels $ of 
class (i) which can be described in terms of n environmental modes prepared into a 
pure state. It is given by 

F = Sr-^S^, (55) 



which follows by imposing the minimal uncertainty condition (1131) to the n-mode 
covariance matrix (l54l) and by using (128]) . Similarly one can verify that given a pure 
n-modes Gaussian state pe and an 5* e Sp{An, R) (|22)) with an invertible subblock 
S2, then the corresponding BGC satisfies condition (|55l) . The above result can be 
strengthened by looking at the solutions for channels of class (ii) of which the channel 
of class (i) are a proper subset. 

To achieve this goal, let us first note that with the choice we made on a^ = cr4„, 
the two matrices a and /3 of Eq. (|45l) are 2n x 2n covariance matrices for two sets of 
independent n bosonic modes satisfying the uncertainty relations ^ with respect to 
the form a2n- In turn, the matrices 5 and S^ of Eq. (l46l) represent cross -correlation 
terms among such sets. After all, the entire covariance matrix 7^ corresponds to a 
pure Gaussian state. 

They key point is now the observation that in Eq. (1431) . the matrix A couples only 
with those rows and columns of the matrices 5 and P which contain elements ^tn-r/2 
or f(^tn-r/2)'- As far as A is concerned, one could indeed replace the element p~^ 
and /(/i^^) of such matrices with zeros. The only reason we keep these elements, 
in the way they are in Eqs. (l45l) and (|46|) . is to render 'Je the covariance matrix of 
a minimal uncertainty state. In other words, the elements of 5 and /5 proportional 
to /i~^ or /(/i^^) are only introduced to purify the corresponding elements of the 
submatrix a, which is in itself hence a covariance matrix of a mixed Gaussian state. 



^Not mentioning the trivial case of Gaussian unitary transformation which does not require any 
environmental mode to construct a unitary dilation. 
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Suppose then that /i of Eq. (l38l) has (say) the first r'/2 eigenvalues equal to 1, 
i.e., 1^1 = 1^2 = ■ ■ ■ = l^r'/2 = 1 while for j E {r'/2 + 1, ■ ■ ■ , r/2} we have that 
Hj E (0, 1). In this case the corresponding sub-matrix of a associated with those 
elements represent a pure Gaussian state, specifically the vacuum state. Accordingly, 
there is no need to add further modes to purify them. Taking this into account, one can 
hence reduce the number of environmental modes ip^^e that allows one to represent <l> 
as in Eq. (|20l) in term of a pure state pe from 2n to 



'pure 



n+{n- r' 12) = 2n - r'/2 , (56) 



i.e. we need the n modes of a plus n — r'/2 additional modes of P to purify those of 
a which are not in a pure state already. An easy way to characterize the parameter r' 
is to observe that, according to Eq. (l37l) . it corresponds to the number of eigenvalues 
having modulus 1 of the matrix of 0T.'0^, i.e., 

r' = 2n- rank[]l2„ - OS'(S')^0^] = 2n - rank[l2n - S'(S')^] 

= 2n - rank[F - S r "^ S^] . (57) 

The explicit expressions for corresponding values of 7^; and S2 are given in Ap- 
pendix IC.ll Here we notice that for r' = r = 2?2 we get £pure = n. This should 
correspond to the channels (|55] ) of class (i) for which one can construct a unitary di- 
lation with pure input states. Indeed, according to Eq. (l57l) . when r' = 2n the matrix 
y — SF^^S^ must be zero, leading to the identity (|55l) . 

Taking into account that r' ^ r = rank[S], a further reduction in the number 
of modes i can be obtained by dropping the requirement of 'Je being a minimal 
uncertainty covariance matrix. Indeed, an alternative unitary representation ([201) of 
$ can be constructed with only 

i = n+{n- r/2) = 2n - r/2 , (58) 

environmental modes (see Appendix IC . 2 1 for the explicit solution). 

The whole analysis can be finally generalized to the BGCs of class (iii), corre- 
sponding to channels that have non invertible matrices Y. We have seen in fact that, 
in this case, the state pe which provides us the unitary dilation of Theorem 1 is con- 
structed by replacing Y with the strictly positive operator Y of Eq. (|49l ). Therefore 
for these channels fpure of Eq. (|56l ) is defined by Eq. (1571 ) with Y replaced by Y, i.e. 

r' = 2n~ rank[F - E Y'^ E^] . (59) 

Taking F := F + (]l2„ — H) with 11 being the projector on Supp[F] this gives, 

r' = 2n-rank[F-SF®^S^] = 2n-rank[F-SrQ^S^] -rank[]l2n-n] 

= fc - rank[r - S r®^ S^] , (60) 

where k = rank[F] = rank[n], where F®^ := nF~^n denotes the Moore-Penrose 
inverse [|26l of Y, and where we have used the fact that Supp[S] C Supp[F]. Re- 
membering then that for channels of class (ii) k = 2n and F®^ = Y^^ these results 
can be summarized as follows: 
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Theorem 2 (Dilations of BGCs involving fewer additional modes) Given $ a BGC 
described by matrices X and Y satisfying the conditions 4771) and characterized by 
the quantities 

r = rank[T] , r' = rank[Y] - rank[Y - SF^^S^] . (61) 

Then it is possible to construct a unitary dilation d20l) of Stine spring form (i.e., involv- 
ing a pure Gaussian state pe) with at most (.pure = 2n — r'/2 environmental modes. 
It is also always possible to construct a unitary dilation f l20l) using i = 2n — r/2 
environmental modes which are prepared in a Gaussian, but not necessarily pure 
state. 

It is worth stressing that, for channel of class (ii) and (iii), the Theorem 2 only 
provides upper bounds for the minimal values of i and ipure- Only in the generic case 
(i) these bounds coincide with the real minima. 

2.3 Minimal noise channels 

In a very analogous fashion to the extremal covariance matrices corresponding to 
pure Gaussian states, one can introduce the concept of a minimal noise channel. In 
this section we review the concept of such minimal noise channels [18] and provide 
criteria to characterize them. Given X,Y E j^2nx2n satisfying the inequality (fT5l) , 
any other Y' = Y + AY with AY ^ will also satisfy such condition, i.e., 

r ^r^2((T2„-xV„x). (62) 

Furthermore, due to the compositions rules (fTSl ). the BGC $' associated with the 
matrices X, Y' can be described as the composition 

$' = ^ o $ , (63) 

between the channel $ associated with the matrices X, Y, and the channel \E' de- 
scribed by the matrices X = !„ and Y = AY. The latter belongs to a special case 
of BGC that includes the so called additive classical noise channels [fTTl |3l HI - see 
Sec. iMlfor details. 

For any X E ]R2"x27i^ ^^^ ^^^ ^^yqu ask how much noise Y it is necessary to add 
in order to obtain a map satisfying the condition (flSl) . This gives rise to the notion 
of minimal noise ifTSl . as the extremal solutions Y of Eq. (fT5l) for a given X. The 
corresponding minimal noise channels are the natural analogue of the Gaussian pure 
state and allows one to represent any other BGC as in Eq. (l63l) with a proper choice 
of the additive classical noise map \1'. 

Let us start considering the case of a generic channel $' of class (i) described by 
matrices X and Y'. According to Theorem 1 it admits unitary dilation with £ = n 
modes described by some covariance matrix 7^ satisfying the condition 

r = S27E4, (64) 

for some proper 2n x 2n real matrix S2. According to Eq. (fTTl) 'Je can be written as 

j'e = je + A, (65) 
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with A ^ and 7^ minimal uncertainty state. Therefore writing Y = S2'JeS2 and 
AY = S2AS2 we can express $' as in (l63l) . where now $ is the BGC associated with 
the minimal noise environmental state 'Je- Most importantly since the decomposi- 
tion (1651 ) is optimal for 7^, the channel $ is an extremal solution of Eq. (fTSl ). We 
stress that by construction $ is still a channel of class (i): in fact it has the same S as 
$', while Y is still strictly positive since 7s > and S2 is invertible - see Eq. (|64l ). 
We can then use the results of Sec. l2.2l to claim that $ must satisfy the equality (|55l) . 
This leads us to establish three equivalent necessary and sufficient conditions for 
minimal noise channels of class (i): 

(66) 
(67) 
(68) 

with r and r' as in Eq. (|6TI) . Since for class (i) we have that r = 2n, the minimal noise 
condition 777,3 simply requires the eigenvalues of the matrix yU of Eq. (|37l ) to be equal 
to unity. Similarly, minimal noise channels in case (ii) and (iii) can be characterized. 

Theorem 3 (Minimal noise condition) A Gaussian bosonic channel characterized 
by the matrices Y and X G ]R2nx2n ^^ ^ minimal noise channel if and only if 



K) 


r = sr-^s^, 


(^2) 


Det[F] = Det[S; 


("^3) 


r = r' , 



Y = SF®^S^, 



(69) 



where, as throughout this work, S = a2n — X a2nX. 



Proof: The complete positivity condition (fT5l) of a generic BGC is a positive semi- 
definite constraint for the symplectic form S, corresponding to the constraint 7 — 
^cr2n ^ in case of covariance matrices of states of n modes. In general, r = rank[S] 
is not maximal, i.e., not equal to 277. When identifying the minimal solutions of the 
inequality (fT?]) . without loss of generality we can look for the minimal solutions of 



Y' - iL' > 0, 



where here 



S' 






/^ 




-/^ 









_ 



(70) 



(71) 



with /i > being diagonal of rank r/2 ( here Y' = OYO^ and S' = OT.O'^ with 
O E 0(277, R) orthogonal). The minimal solutions of inequality (iTOl) are then given 
by Y' = SS'^ © 0, where S* is a r x r matrix satisfying 



S 



/i 
-/i 



s' 





-/i 







(72) 



so a symplectic matrix with respect to the modified symplectic form, so an element 

of {Af G Gl{r,R) : M = (/iV2 © ^1/2)5(^-1/2 ^ ^-1/2)^ ^ ^ Sp{r,R)}. From 
this, it follows that the minimal solutions of (iTOl) are exactly given by the solutions 
of Y' = S'(F')®^ (S')^, from which the statement of the theorem follows. ■ 
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2.4 Additive classical noise channel 

In this subsection we focus on the maps \& which enter in the decomposition (1631) . 
They are characterized by having X = l2n and F ^ 0. Note that with this choice 
the condition (flSl) is trivially satisfied. This is the classical noise channel that has 
frequently been considered in the literature (for a review, e.g., see Ref. JH). For 
completeness of the presentation, we briefly discuss this class of multi-mode BGCs. 
If the matrix Y is strictly positive, the channel ^ is the multi-mode generalization 
of the single mode additive classical noise channel flTlElllll. In the language of Ref. 
ifTTl . these are the maps which have a canonical form B2. Indeed, one can show 
that these maps are the (Gaussian) unitary equivalent to a collection of n single mode 
additive classical noise maps. To see this, let us apply symplectic transformations (5*1 
and S2) before and after the channel ^. Following Eq. (fT9l ) this leads to {1„, Y} 1— »• 
{S'iS'2, 5^1^52}. Now, since F > 0, according to Williamson's theorem [|23l , we 
can find a 6*2 G Sp{2n, R) such that S2YS2 is diagonal diag(Ai, ■ ■ ■ , A„, Ai, • • • , An) 
with Ai > 0. We can then take Si = 82^ to have S'iS'2 = t2n- For F ^ but not 
Y > 0, the maps \1/ that enter the decomposition Eq. (|63l) however include also 
channels which are not unitarily equivalent to a collection of B2 maps. An explicit 
example of this situation is constructed in Appendix |El 

2.5 Canonical form for generic channels 

Analogously to Refs. [fTTl [T3l [T4ll . any BGC $ described by the transformation 
Eq. (flTl) can be simplified through unitarily equivalence by applying unitary canoni- 
cal transformations before and after the action of the channel which induces transfor- 
mations of the form (fT9l) . Specifically, given a n-mode Gaussian channel <l> described 
by matrix X and Y we can transform it into a new n-mode Gaussian channel $c de- 
scribed by the matrices 

X, = S1XS2 , Fc = SJYS2 , (73) 

with Si,2 G Sp{2n, R). As already discussed in the introductory sections, from an 
information theoretical perspective $ and $c are equivalent in the sense that, for in- 
stance, their unconstrained quantum capacities coincide. We can then simplify the 
analysis of the n-mode Gaussian channels by properly choosing Si and S'2 to induce a 
parametrization of the interaction part (i.e., X) of the evolution. The resulting canon- 
ical form follows from the generalization of the Williamson theorem [|23l presented 
in Ref. [|20l . According to this result, for every non-singular matrix X E Gl{2n, R), 
there exist matrices S'i_2 G Sp{2n, R) such that 



X(, — blX S2 



In 

J^ 



(74) 



where J^ is an x n block-diagonal matrix in the real Jordan form f26\ . This can be 
developed a little further by constructing a canonical decomposition for the symplec- 
tic matrix S associated with the unitary dilation (|20|) of the channel. 

For the sake of simplicity in the following we will focus on the case of generic 
quantum channels $ which have non-singular X E Gl{2n, R) and belong to the class 
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(i) of Theorem 1 (i.e., which have r = rankfS] = 2n). Under these conditions X 
must admit a canonical decomposition of the form (fT?!) in which all the eigenvalues 
of J are different from 1. In fact one has 

S = 0"2n ~ X a2nX = 5*2 [o'2n ~ X^ a2nXc\ 5*2 = 5*2 Sc 5*2 , (75) 

with Sc being the skew-symmetric matrix associated with the channel $£, i-e-, 



E. : = 



tn- J 

r - K 



(76) 



Since rank[Sc] = rank[S] = 2n, it follows that J cannot have eigenvalues equal to 
1. Similarly, it is not difficult to see that if X has a canonical form (1741) with all the 
eigenvalues of J being different from 1, then $ and $c are of class (i). However, a 
special case in which X = ]l2n is investigated in details in Appendix |El 

Consider then a unitary dilation (|20l ) of the channel <l>c constructed with a not 
necessarily pure Gaussian state pe of i = n environmental modes. According to 
the above considerations, such a dilation always exists. Let S E Sp{An, H) be the 
An X An real symplectic transformation (l22l) associated with the corresponding unitary 
U. Assuming si = Xj, an explicit expression for this dilation can be obtained by 
writing 



S4 



In 






J' 











(77) 



where, for j = 2, 3, Fj, Gj are n x n real matrices. Imposing Eqs. (1231) . one obtains 
the following relations 



r + F2GI 
Gl + F2J' 



IT 



0, 



J'^ + F,Gl 
Gl + F3 J^ = 



1, 



(78) 



whose solution gives an S' G Sp{An^ R) of the form 



S 



In 






J 



-G^J-^ 





G-^J{J 



't. 



r)G-^ 






G 



fln 





G-^JG 



(79) 



with G being an arbitrary matrix G E Gl{n, R). As a consequence of this fact, and 
because the eigenvalues of J are assumed to be different from 1, S2, S3 and S4 are 
also non-singular. This is important because it shows that in choosing S as in the 
canonical form (|79l ) we are not restricting generality: The value of S2 can always be 
absorbed into the definition of the covariance matrix 'Je of Pe by writing 



lE 



-1 



Yrs: 



(80) 



(see also Appendix ID]). Taking this into account, we can conclude that Eq. (1791 ) 
provides an explicit demonstration of Lemma 1 for channels of class (i) with non- 
singular X. 
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Since <l>c is fully determined by Xc and Yc, the above expressions show that the 
action of $c on the input state does not depend on the choice of G. As a matter of 
fact, the latter can be seen as a Gaussian unitary operation Uq characterized by the 
n-modes symplectic transformation Sp{2n, R), 



A 



G 



\ G^ 








G~' \ 



(81) 



applied to final state of the environment after the interaction with the input, i.e., 
^G = Ug^Uq, where $ is the weak complementary map for G = !„, and <i>G is the 
weak complementary map in presence of G 7^ !„ - see the next section for details. 
Since the relevant properties of a channel (e.g., weak degradability |fT6l[T7l ) do not 
depend on local unitary transformations to the input/output states, without loss of 
generality, we can consider G = ~J and the canonical form for S E Sp{An, R) 
assumes the following simple expression 



S 





J 









In- J 



- J- 







-J 



In 





J 



(82) 



The possibility of constructing different, but unitarily equivalent, canonical forms for 
S is discussed in Appendix iDl 



3 Weak degradability 

Among other properties, the unitary dilations introduced in Section |2] are useful to 
define complementary or weak complementary channels of a given BGC $. These 
are defined as the CPT map $ which describes the evolution of the environment under 
the influence of the physical operation describing the channel [i16li17iL i.e.. 



^{p):=1xs[U{p®Pe)U^] 



(83) 



where p, pe and (j are defined as in Eq. (|20l ), but the partial trace is now taken over 
the system modes. 

Specifically, if the state pe we employed in constructing the unitary dilation of 
Eq. (|20l) is pure, then the map $ is said to be the complementary of $ and, up to 
partial isometry, it is unique Il27ll28ll29l[30ll . Otherwise it is called weak complemen- 
tary liT6l[T7l . Since in Eq. (|20|) the state pE is Gaussian and f/ is a unitary Gaussian 
transformation, one can verify that $ is also BGCcl. Expressing the Gaussian unitary 
transformation U in terms of its symplectic matrix S of Eq. (|22)) the action of <l> is 
fully characterized by the following mapping of the covariance matrices 7 of p, i.e.. 



$: 



7 



■537^ + Si-^ESI 



(84) 



^In general however, it will not map the n input modes into n output modes. Instead it will 
transform them into £ modes, with I being the number of modes assumed in the unitary dilation ( |20l ). 
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which is counterpart of the transformations (fT6l) and (l24l) that characterize $. The 
channel $ is then described by the matrices X = s^ and Y = s^'JesJ which, ac- 
cording to the symplectic properties (|23l) . satisfy the condition 



y ^ iS with S := af^ - X^(T2nX . (85) 

The relations between $ and its weak complementary $ contain useful informa- 
tion about the channel $ itself. In particular we say that the channel $ is weakly 
degradable (WD) while $ is anti-degradable (AD), if there exists a CPT map T 
which, for all inputs p, allows one to recover <i>(p) by acting on the output state $(/3), 
i.e. 

T o $ = 1) . (86) 

Similarly, one says that $ is AD and $ is WD if there exists a CPT map T such that 

T o $ = $ . (87) 

Weak degradability lfT6l[T7ll is a property of quantum channels $ generalizing the 
degradability property introduced in Ref. [|27l . The relevance of weak-degradability 
analysis stems from the fact that it allows one to simplify the quantum capacity sce- 
nario. Indeed, it is known that AD channels have zero quantum capacity [f6l[T7l. 
while WD channels with pE pure are degradable and thus admits a single letter ex- 
pression for this quantity ETl . A complete weak-degradability analysis of single 
mode bosonic Gaussian channels has been provided in Ref. lfT6l[T7l . Here we gener- 
alize some of these results to n > 1. 

3.1 A criterion for weak degradability 

In this section we review a general criterion for degradability of BGCs which was 
introduced in Ref. [[TSl . adapting it to include also weak degradability. Before en- 
tering the details of our derivation, however, it is worth noticing that generic multi- 
mode Gaussian channels are neither WD nor AD. Consider in fact a WD single-mode 
Gaussian channel $ having no zero quantum capacity Q > (e.g., a beam-splitter 
channel with transmissivity > 1/2). Define then the two mode channel $ C?) $ with 
$ being its weak complementary defined in fll6[[T7l . This is Gaussian since both $ 
and $ are Gaussian. The claim is that $ (g) $ is neither WD nor AD. Indeed, its weak 
complementary can be identified with the map $ C?) $. Consequently, since $ (g) $ 
and $ ® $ differ by a permutation, they must have the same quantum capacity Q' . 
Therefore if one of the two is WD than both of them must also be AD. In this case Q' 
should be zero which is clearly not possible given that Q' ^ Q. In fact, one can use 
$ ® $ to reliably transfer quantum information by encoding it into the inputs of $. 
In this respect the possibility of classifying (almost) all single-mode Gaussian maps 
in terms of weak degradability property turns to be rather a remarkable property. We 
now turn to investigating the weak degradability properties of multi-mode bosonic 
Gaussian channels deriving a criterion that will be applied in Sec. |4.1| for studying in 
details the two-mode channel case. 
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Consider a n-mode bosonic Gaussian channel $ characterized by the unitary di- 
lation (|20l) and its weak complementary $ (l83l) . Let {X, Y}, {X, Y} be the matrices 
which define such channels. For the sake of simplicity we will assume X and X to 
be non-singular, with X, X G Gl{2n, R). Examples of such maps are for instance 
the channels of class (i) with X non-singular described in Sec. 12.51 Adopting in fact 
the canonical form (f82l ) for S we have that 



X 











X 



K 









(88) 



with all the eigenvalues of J being different from 1. 

Suppose now that $ is weakly degradable with T being the connecting CPT map 
which satisfies the weak degradability condition (|86l) . As in Refs. [ 161 [TTll we will 
focus on the case in which T is BGC and described by matrices {Xr, Ir}- Under 
these hypothesis the identity (f86l ) can be simplified by using the composition rules 
for BGCs given in Eq. (fTSi) . Accordingly, one must have 



Xr 
Yr 



x-^x, 

Y - X^YXt . 



(89) 



These definitions must be compatible with the requirement that T should be a CPT 
map which transforms the n system modes into the i environmental modes, i.e.. 



Yt ^ % (ai - X^(T2„Xr) 



(90) 



Combining the expressions above, one finds the following weak-degradability con- 
dition for n-mode bosonic Gaussian channels hlSl . i.e. 



Y - X^X-^(F + io2n)X-^X + io'^f ^ , 



(91) 



In order to obtain the anti-degradability condition (|87l) . it is sufficient to swap {X, F} 
with {X, F} and the system commutation matrix oin with a^, in Eq. (|9T1) . i.e., 

Y - X^X-^(r + io^f)X-^X + ia2n ^ . (92) 

Equations (|9T1) and (|92l) are strictly related. Indeed since 

Y - X^X-^(f + ia^f)X-^X + ia2n (93) 

= -X^X"^ (y - X^X-^{Y + i(j2n)X~^X + ia^^ X~^X , 

equation (|92l ) corresponds to reverse the sign of the inequality (|9T| ). i.e. 

F - X^X-^(F + ia2n)X~^X + iaf^ ^ . (94) 

Hence to determine if $ is a weakly degradable or anti-degradable channel, it is then 
sufficient to study the positivity of the Hermitian matrix 



W :=? - X^X-^{Y + ia2n)X-^X + iaf^ . 



(95) 
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In the case in which £ = n this can be simplified by reminding that an Hermitian 
2n X 2n matrix W partitioned as 



W 



Wi W2 



with Wi being n x n matrices is semi-positive definite if and only if 
VTi ^ and Ws - W^W{^W2 ^ , 



(96) 



(97) 



the right hand side being the Schur complement of W (see, e.g., page 472 in Ref. 
Il26l ). Using this result and the canonical form (|82l) . Eq. (|9T| ) can be written as in 
Eq. dlT]) with 



Wl 

W2 
W^ 



l{J-'^ 



j-^)-'YAtn-J-')-' 



21, 



Y2{J 



-T 



Yi 



(98) 



J-^)-'Y, 



Y; - (j-i - tn)Y^{r 



and 



Y 



Yi Y2 
Y^ Y, 



(99) 



For the anti-degradability condition (l92l) simply replace [^] with [^] in Eq. (|97 



4 Two-mode bosonic Gaussian channels 

Here we consider a particular case of n-mode bosonic Gaussian channel analysis 
above, namely, the case of n = 2. This is by no means such a special case as one 
might at first be tempted to think since any ra-mode channel can always be reduced 
to single-mode and two-mode parts [|20ll . For two-mode channels the interaction 
part and the noise term of a generic two-mode bosonic Gaussian channel, X and 
Y , respectively, are 4x4 real matrices. Particularly, we will focus on two-mode 
channels $ which have non-singular X and belong to the class (i) of Theorem 1 (i.e., 
which have r = rank[S] = 4), like in Sec. 12.51 These maps can be grasped in terms 
of a unitary dilation of the form (f82l) coupling the two system bosonic modes with 
two additional (environmental) modes, where J is a 2 x 2 real Jordan block. In order 
to characterize this large class of two-mode BGCs, one has to examine only three 
possible forms of J: 

• Class A: This corresponds to taking a diagonalizable Jordan block, that is. 



J:= Jn 



a 
h 



(100) 



where a and h are real nonzero numbers. It represents the trivial case of a two- 
mode bosonic Gaussian channel, whose interaction term does not couple the 
two modes. Actually, we call it of class Aiiia ^h and of class A2 otherwise. 
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Class B: This is to take J as a non-diagonalizable matrix with a nonzero real 
eigenvalue a with double algebraic multiplicity (but with geometric multiplic- 
ity equal to one), i.e. 



J:= Ji 



a 1 
a 



(101) 



In this case the Jordan block is called defective [l26l|. Here, a noisy interaction 
between the bosonic system and the environment, coupling the two system 
modes, is switched on. 



Class C: Here the real Jordan block J has complex eigenvalues, i.e. 

J:=J2-- 



a h 
—h a 



(102) 



with & 7^ 0; the eigenvalues of J are a ± ih. Again, the two system modes are 
coupled by the noisy interaction with the environment through the presence of 
the element h. 

In order to explicit the form of F = S27£;s|^, with S2 being defined as in Eq. (|82l ). 
we consider a generic two-mode covariance matrix in the so-called standard form 
[|32l for the environmental initial covariance matrix 7^, i.e. 



lE 



where 



Ti, 



Ti 








^2. 


X 


Z^^ 


Z-,+ 


y 



2„,2 



(103) 



(104) 



r- 



and X, y, z+- are real number satisfying x + y ^ 0, xy — z_ ^1 and x y 
x"^ + (z--2+ — 1)^ — xy{z'^ + z'^) ^ because of the uncertainty principle. More 
generally, one can apply a generic two-mode (symplectic) squeezing operator V{e) 
to the environmental input state, i.e.. 



7'e 



V{e)^EV{e)' 



(105) 



V{e) 







where 

„ _ c + hs —qs 

R ' —qs c — hs 

and c = cosh(2r), s = sinh(2r), h = cos(20), q = sin(20) and e = re 
squeezing parameter [|32l . Finally, it is interesting to study how the canonical forms 
of two-mode BGCs compose under the product. A simple calculation shows that the 
following rules apply 



(106) 
2^"^ being the 



A 



B 



C 



A 
B 
C 



A 

Ai/B 
Ai/B/C 



Ai/B 

A2/B 
AJB/C 



Ai/B/C 

Ai/B/C 

A/C 
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In this table, for instance, the element on row 1 and column 1 represents the 
class (i.e.. A) associated to the composition of two channels of the same class A. 
Note that the canonical form of the products with a "coupled" channel (i.e., with B 
or C) is often not uniquely defined. For instance, composing two class B channels 
characterized by the matrices 



(J 



li 



ai 1 
a,: 



(107) 



with 2 = 1,2, will give us either a class A2 channel (if ai + 02 = 0) or a class B 
channel (if ai + 02 7^ 0). Composition rules analogous to those reported above have 
been analyzed in details for the one-mode case in Ref. [[TtI . In the following we will 
study the weak-degradability properties of these three classes of two-mode Gaussian 
channels. 

4.1 Weak-degradability properties 

The weak-degradability conditions in Eqs. (1971 ) become 

Ti - (]l2 - r^)Ti{t2 - J'^) ^ (108) 

and 

jr2J^-(ii2- J)r2(ii2- J^) (109) 

-( J-i - 2112) [Fi - (I2 - J~^)Ti{t2 - J-^)] "' {J-^ - 2U) ^ . 

In the same way, the anti-degradability is obtained when both these quantities are 
non-positive. As concerns the environmental initial state of the unitary dilation, one 
can consider a generic two-mode state as in Eq. (|105l) . On one hand, we find that, if 
[J, R\ = 0, this two-mode squeezing transformation V{e) can be simply "absorbed" 
in local symplectic operations to the output states and then it does not affect the 
weak-degradability properties. On the other hand, if [J, R] ^ 0, we find numerically 
that the introduction of correlations between the two environmental modes contrasts 
with the presence of (anti-) weak-degradability features. Therefore, one can consider 
the particular case in which the environment is initially in a state with a symmetric 
covariance matrix 7^ as in Eq. (11031) with x = y = 2N + 1 and z^ = z+ = 
where A^ ^ 0. In this case je = (2A^ + 1)12 corresponds to a thermal state of two 
uncoupled environmental modes with the same photon average number A^ and it is 
possible to see the results above easily through analytical details. In fact, we study 
analytically the positivity condition in Eq. (|9T| ) in the three possible forms of the real 
Jordan block Jj. 

In the uncoupled case Jq as in Eq. (|100l) . substituting in Eq. (|9T1) , we find that 
these two-mode bosonic Gaussian channels are WD if a,b ^ 1/2 and AD for a,b < 
1/2 (any A^ ^ 0). In other words, in the case of two uncoupled modes, the weak- 
degradability properties can be derived from the results for one-mode bosonic Gaus- 
sian channels: tensoring two WD (AD) one-mode Gaussian channels with WD (AD) 
one-mode Gaussian channels yield two-mode Gaussian channels which are WD (AD). 
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In the case of defective J, i.e., Ji as in Eq. (|101l) . corresponding to noisy inter- 
action coupling the two system modes, substituting in Eq. (|9T| ). we find that, on one 
hand, these two-mode bosonic Gaussian channels are WD if a > 1 and 



AT ^ A^i := - 
^ 2 



\2a 



1 

2 ^/a{a 



(110) 



On the other hand, it is AD if a < and N ^ Ni (see Fig. [21) . Note that the defective 
Jordan blocks are not usually stable with respect to perturbations f20\. Indeed, we 
find numerically that, applying proper two-mode squeezing transformations to the 
environmental input, these weak-degradability conditions reduce to the decoupled 
case ones. In Fig. [T]we consider, for simplicity, a symmetric environmental initial 
state 7^ as in Eq. (11051) with x = y, z^ = and e = r, and we plot the relation 
between x, z+ and the minimum value of r such that J := Ji reduces to J := Jq 
corresponding to the decoupled case. One realizes that a squeezing parameter r close 
to 1 is enough to decouple the two modes representing the system, carrying quantum 
information. Moreover, let us point out that this squeezing threshold (r) increases 
slightly with the presence of correlations (z^) while decreases when increasing the 
level of noise (x) in the initial environmental state 7^. 




Figure 1 : Relation between the parameters x, z^ and the minimum value of r in the 
initial environmental state such that the two-mode channel with X = I2 © Ji reduces 
to the decoupled case X' = 1 2 © Jo with the same interaction parameter a for the 
two system modes. 



Finally, in the case of real Jordan block with complex eigenvalues, i.e., J2 as 
in Eq. (|102k the corresponding two-mode bosonic Gaussian channels are WD if 
a > 1/2 and 



N> No 



462 



2a)2 



1/2 



(111) 



while they are AD if a < 1/2 and N ^ N2 (see Fig. |2]). In both of these cases (real 
and complex eigenvalues), in which the interaction term couples the two bosonic 
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Figure 2: In continuous line we report A^^i as function of a in the case of Ji. For 
N ^ Ni the map is WD if a > 1 and AD if a < 0. In dashed line we plot N2 as 
function of b when a = in the case of J2. For N > N2 the channel is WD (AD) if 

a > 1/2 (a < 1/2). 



modes, there is the (apparently) counter- intuitive fact that above a certain environ- 
mental noise threshold the weak-degradability features appear, while for one-mode 
bosonic Gaussian channels they do not depend on the initial state of the environ- 
ment. Actually, one would expect at most that, when the level of the environmental 
noise increases, the coherence progressively decreases until to be destroyed. It would 
mean that it becomes more and more difficult to recover the environment (system) 
output from the system (environment) output after the noisy evolution. However, 
the things go the other way around when multi-mode bosonic Gaussian channels are 
considered. 

4.2 Channels with zero quantum capacity 

Analogously to Ref. [|T71 where the one-mode case is investigated, one can enlarge 
(other than the AD maps) the class of two-mode BGCs with Q = 0, composing a 
generic channel with an AD one. First of all, consider a channel $ as in Section 
12. 3[ but being AD (not necessarily minimal noise), then the maps $', defined in 
Eq. (|63l) . have zero quantum capacity, i.e., they cannot be used to transfer quantum 
information. For instance, one can choose 7E = {2Nc + 1) !„, i.e., the environmental 
initial state of the map $ is a multi-mode thermal state with N^ being the average 
photon number for each mode, such that $ is AD or simply with zero capacity; 
therefore, for any 7^ ^ 75; = (2iVc + 1)]1„, as in Eq. (l65l) . the map $' of Eq. (|63l) 
has Q = 0. Particularly for n = 2, using these observations and choosing Nc equal 
to either Ni (and a < 0) or N2 (and a < 1/2) as in Eqs. (IllOl ) and (lllll) . one 
obtains that for X = U 2 © Ji,2 and Y' = S27^s^ [with S2 as in Eq. (|82l )1 the resulting 
channel $' has always zero capacity. In this way, one extends considerably the set of 
two-modes maps with zero capacity, other than the very particular cases of two-mode 
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environmental thermal states studied above and shown in Fig. [2l For instance, two- 
mode squeezing can be applied to the thermal state 7^ including not only states with 
N > Nc but also with not trivial two-mode correlations such that 7^ ^ (2iVc + l)]l2- 
Therefore, just considering this last simple inequality one includes so a larger set of 
maps that have zero quantum capacity. 

Moreover, we observe that, according to composition rules above, the combina- 
tion $ = $// o $j of two channels $/ and $// of class A2 and C, respectively, with 
Jordan blocks Jj as in Eq. (1 1001) with aj = bj and Jjj as in Eq. (11021) with ajj and 
bn 7^ 0, gives J = ajJn which is in the class C. Now, since we have Ni ^ 0, 
N2 ^ and assuming a/ ^ 1/2, the channel $/ is AD and the resulting channel $ 
must have Q = 0. Varying the parameters but keeping the product ajajj = a and 
o-ibii = b fixed, the parameter A^ can assume any value satisfying the inequality 



-^i 



5(l-4a + 8a2 + 862)y/2^ 
62 + (g-1)2 



(112) 



Note that a/ has been chosen equal to 1/2 and $/ corresponds to two uncoupled 
beam-splitter maps with transmissivity 1/2. We can therefore conclude that all chan- 
nels of the form C with N as in Eq. (11121) have zero quantum capacity - see Fig.[3l 

Consider now the composition $ = $// o $/ of two channels $/ and $// of class 
C and A2 (i.e., in the opposite order with respect to above), respectively, with Jordan 
blocks Jj as in Eq. (|102l) with aj and bj ^ and Jji as in Eq. (|100l) with aji = bn, 
giving J = aiiJi which is in the class C . As before, since we have A^i ^ 0, A^2 ^ 
and assuming again an ^ 1/2, the channel $2 is AD and the resulting channel has 
(5 = 0. Varying the parameters but keeping the product a/O// = a and b^au = b 
fixed, the parameter N can assume any value satisfying the inequality 



N > i 

4 



;i + 4a2 + 4&2) (1 - 4a + 8a^ + 8b^) \ '/' _ 
4(62-f(a-l)2)(a2 + 62) 



(113) 



where again an is chosen equal to 1/2. Again we can conclude that all class C 
channels with A^ as in Eq. (11131) have zero quantum capacity. However, notice that the 
constraint in Eq. (|113l) is an improvement with respect to the constraint of Eq. (|112l) 
- see Fig. |3l 

5 Conclusions 

In this work, we have presented a complete analysis of generic multi-mode Gaus- 
sian channels by proving a unitary dilation theorem and by finding their canonical 
form. This is a simple form that can be achieved for any Gaussian quantum channel, 
as a convenient starting point for various considerations. For instance, it allows us 
to simplify the analysis of the weak-degradability properties of multi-mode bosonic 
Gaussian channels. Minimal output entropies, or quantum and classical information 
capacities and other difficult questions might be tackled using the canonical form of 
multi-mode Gaussian channels shown in this paper. Here, we investigated in details 
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Figure 3: The continuous line depicts plot N2 as in Eq. (|111|) versus h, with a = 1 
in J2 of Eq. (IT02l) . For A^ > A/^2 the channel is WD (AD) if a > 1/2 (a < 1/2). 
The dashed line refers to the bound in Eq. (11121) . while the dashed-dot line to the 
one in Eq. (11131 ); above these bounds the class C map is WD but with Q = 0. Note 
that Eq. (11131 ) is an improvement with respect to the constraint of Eq. (|112l) . Similar 
bounds can be obtained in the case a < 1/2, enlarging the group of AD maps with 
other channels with Q = 0. 



the two-mode scenario that is relevant since any n-mode channel can always be re- 
duced to single-mode and two-mode parts ||20l . Furthermore, the results of this paper 
could play a basic role in characterizing the efficiency of continuous -variables quan- 
tum information processing, quantum communication and quantum key distribution 
protocols. 

6 Acknowledgements 

FC. and V.G. thank the Quantum Information research program of Centro di Ricerca 
Matematica Ennio De Giorgi of Scuola Normale Superiore for financial support. I.E. 
acknowledges the EPSRC, the QIP-IRC, the EU (QAP, COMPAS), Microsoft Re- 
search, and the EURYI Award Scheme for financial support. A.S.H. acknowledges 
support of RFBR grant 06-01 -00 164-a and the RAS program "Modem problems of 
theoretical mathematics". 



A Proof of Lemma 1 

Note that it does not restrict generality to take cr^ = a2i, as this can always be ac- 
companied by an appropriate similarity transform. Our problem at hand of extending 
a symplectic form is then equivalent to the following problem: Suppose we are given 
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column vectors ei, ■ ■ ■ , e„ and /i, ■ ■ ■ , fn from 11^*^"+^) that satisfy 

eja2{n+e)ek = 0, (114) 

fj (^2(n+i)fk = 0, (115) 

ej(T2{n+£)fk = Sj,k, (116) 

for j,k = 1, • • • , n. The procedure continues by identifying vectors e„+i and /„_|_i 

such that el^^a2(n+i)fn+i = 1 and 






0-2(n+^)W^ = /J+lO-2(n+£)W 



for all 
Now define 



w e Wn := span(ei, ■ ■ ■ , e„, /i 







1 Jn, 



(117) 

(118) 
(119) 



W;^ = {w: w'^a2{n+e)V = 0\/v e Wn}. 
It is now not difficuh to see that WnOW^ = {0} and R2{n+^) ^yy^Q^yy^. Suppose 
that the vector v has v'^a2(n+e)<^j ='■ C(j and v^a2{n+e.)fj ='■ Pj ior j = 1, ■ ■ ■ ,n. 
Then 

n 1 r n 

+ t; + ^(a,/,-/5,e,) , (120) 

i=i 

where the first term is element of Wn and the second of W^. Following a symplectic 
Gram-Schmidt procedure, the symplectic basis can hence be completed, which is 
equivalent to extending the matrices si and S2 to a symplectic 

Sl S2 






S 



S3 S4 



G5p(2(n + £),R). 



(121) 



B Derivation of Eq. (1351) 



Here we show that Eq. (|35l) admits solution for s'2 as in Eq. (|4T]) . In fact, assuming 
0"!^ = (^2n © cr2n with cr2„ as in Eq. ©, one has 



4^fn(4)^ 



s' = [ i^-i I O^A ] 



0"2n 












0"2n _ 



S' 



= ir- Va^fT-^ + O'^A G2n A^'O - S' 

= K-^ [KYlK^ + B) K~'^ + O^A a2„ J^ O - Yl 

= K-^BK~^ + 0^Aa2nA^0 

= 0{M'/^BM'/^ + Aa2nA^)0^, (122) 

where we used Eq. (l40l) to write (72„ = KH'K^ + B, with i? being the 2n x 2n 
matrix 



E 





3 











fln-r/2 








( 


) 





-lln-r/2 



(123) 



The identity (l35l) finally follows by noticing that the last term in Eq. (11221) cancels 
since M^^^B = BM^/"^ = 5 and A a2n A^ = -B. 
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C Properties of the environmental states 



In this appendix we first give an explicit derivation of Eq. (I43I) . Then we analyze in 
details the property of the state pE associated with the co variance matrix 7^; defined 
be the Eqs. ^B and ^B- Replacing Eq. dH]) into Eq. ^1, we get 



a 



2n 



S2IE' 



J\T 



a 


6 ' 






[s'^ 


(^\ 



[ K~' I O^A ] 

R-^ a K-^ + O^ A6^ R-^ + K'U A^ O + O^ A (3 A^ O 

O^ (M^/^ a M^''^ + A5'^ M^'^ + M^'"^ 5A^ + A[5A^)0, 



which leads to 



M-^ = a + M-^/'^A 5'^ + 5 A^ M-^'^ + M'^/^ A (3 A^M-^'"^ , 



(124) 



and hence to Eq. (gS]) by the fact M'^/'^A = A^ M'^''^ = A = A^ . Such an 
equation admits the solution given in Eqs. (l45l) and (|46l) . Explicitly this corresponds 
to the An x An covariance matrix 7^; of the form 



r /i"^ 











/(/^-') 








a 





/(a) 





/i-^ 





/(/^-^) 











it 





/(a) 





/(/i-^) 





^^-' 





n 





/(a) 





a 




f{^^-') 











^^-' 








/(a) 





ei . 



where for easy of notation 1 := tn-r/2- By looking at the structure of this covariance 
matrix, one realizes that it is composed by two independent sets formed by r and 
2n — r modes, respectively. The first set describes r/2 thermal states characterized 
by the matrices yU~^ which have been purified adding further r/2 modes. The second 
set instead describes a collection of 2(r2 — r/2) = 2n — r modes prepared in a 
pure state formed hy n — r/2 independent pairs of modes which are entangled. By 
reorganizing its rows and columns this can be cast into the simpler form 



lE 



K' 



fm 



^,-l^ 



fm- 



r' 







a 



2n—r 



fiih^-r) 



/(a 



2n-r) 



ei 



2n—r 



}r 
} r 

} 2r2 — r 
} 2n — r - 



(125) 



where we used fx to indicate the r x r matrix ji = jj,® H- 



C.l Solution for ^pure = 2n — r' /2 environmental modes 

Defining r' as in Eq. (|57] ) we choose the environmental commutation matrix to be 
o"^ = cr2„ © <J2n^r' with 02ri and cr2n-r' as in Eq. ([T]). A unitary dilation with £pure = 
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2n—r'/2 environmental modes in a pure state is obtained by having S2 = Y^^'^s^ with 
s'2 as in Eq. (|4T1) . In this case, however, A is a rectangular matrix 2n x 2(n — r'/2) 
of the form 



A 
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( 













an-r/2 
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ln.-r/2 





]r'/2 

} (r - r')/2 

} n — r/2 

}r72 

} (r - r')/2 

} n — r/2. 



(126) 



Similarly, the co variance matrix ^e can be still expressed as in Eq. (1441) . In this case, 
yet, a is a 2n X 2n matrix of block form 
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lr'/2 













^^-o' 














an-r/2 
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^^o' 











^ln-r/2 . 



}r'/2 

}(r-r')/2 
}n — r/2 
}r'l2 

}(r-r')/2 
}n — r/2, 



(127) 



where ^ = 5/4 and yUo is the (r — r')/2 x (r — r')/2 diagonal matrix formed by the 
elements of /i which are strictly smaller than 1. (3 is the (2n — r') x (2n — r') matrix 



/9 



/^o 







.^]ln.-r/2 







/^o 



} (r - r')/2 
} ra — r/2 
} (r - r')/2 
^an-r/2 J }ri-r/2, 



(128) 



and 
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f{itn^r/2) 
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f{iK-r/2) 





}r72 

}(r-r')/2 
jn — r/2 
}r'/2 
}(r-r')/2 
}n — r/2, 



(129) 



with / as in Eq. (|46l ). 

By looking at the structure of this covariance matrix, one realizes that it is com- 
posed by three independent pieces. The first one describes a collection of r'/2 vac- 
uum states. The second one, in turn, describes (r — r')/2 thermal states characterized 
by the matrices ji^^ which have been purified by adding further (r — r')/2 modes. 
The third one, finally, reflects a collection of 2(ri — r/2) = 2n — r modes prepared 
in a pure state formed by n — r/2 independent pairs of modes which are entangled. 
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C.2 Solution for (, = 2n — r/2 not necessarily pure environmental 
modes 

In this subsection, we present the alternative derivation of a dilation that does not 
necessarily involve an environment prepared in a pure state. Choosing the commu- 
tation matrix a^ = cr2„ © (T2n-r with a2n and o-2n-r as in Eq. ([U), the matrix s'2 can 
be still expressed as in Eq. (|4TI) . In this case, however, A is a rectangular matrix 

2n X (2n — r) of the form 
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ln-r/2 








ln-r/2 



}r/2 

} n — r/2 

}r/2 

} ra — r/2, 



(130) 



Similarly, 7^; has the block form (l44l) . where a is still the 2n x 2n matrix of Eq. (|45] ), 
while /? and 5 are, respectively, the following (2ri — r) x (2n — r) and 2n x (2n — r) 
real matrices: 
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^ln-r/2 











} n — r/2 
^ln-r/2 J }n-r/2, 



(131) 
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}r/2 

}n — r/2, 



with E, and / as in Eq. (|46l ). That is. 



7s 



r /i-^ 











1 





ei 





/(a) 





/i-^ 














ei 


/(ei) 














/(ei) 


a 








/(ei) 











a J 
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}r/2 
} n — r/2 
}r/2 
jn — r/2 
}n — r/2 
}n — r/2 



(133) 



with 1 = Un-r/2- This covariance matrix now consists of two independent parts: The 
first one describes a collection of r/2 thermal states described by the matrices fi^^. 
The second instead reflects a collection of 2(ri — r/2) = 2ri — r modes prepared in a 
pure state formed by n — r/2 independent couples of modes which are entangled. The 
covariance matrix given in Theorem 1 can be recovered from the one given above by 
adding r modes to purify the thermal states /i~^. 



D Equivalent unitary dilations 



Let 



S 



Sl S2 
S3 S4 



(134) 
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and 7s define a unitary dilation for a bo sonic Gaussian channel $ characterized by 
matrices X and Y. Then a full class of unitary dilations 



S' 






(135) 



can be obtained by taking 7^ = V'^eV'^ and 

S[ = Si , 4 = S2V 4 = PFsg , 4 = iyS4l/ , (136) 

with V G Sp{2i, R) and 1^ G Sp{2n, R) being symplectic transformations of i and 
n modes respectively. With this choice in fact 7^ is still a co variance matrix while 
the conditions (l23l) and ((24)) are automatically satisfied. From a physical point of 
view the symplectic transformations V and W correspond to unitary local operations 
applied to the environmental input and output states, respectively, by virtue of the 
metaplectic representation. Consequently, the weak complementary channels $ and 
$' associated with these two representations are unitarily equivalent and the weak- 
degradability properties one can determine for $ will be the same when studied for 

Conversely, let us suppose to have two unitary dilations of $, realized with i = n 
environmental modes and characterized by the symplectic matrices S and S' as in 
Eq. (11341) and (|135l) . respectively, with Sj and s'^ being 2n x 2n square matrices. Then 
it is possible to show that they must be related as in Eq. (11361) under the hypothesis 
that S2 and S3 are non-singular. First of all, since Eq. (I24l) must be satisfied for all the 
input covariance matrices 7, we have si = X"^ = s[. Define then V = Sg^^Sg and 
W = 53^3^^. By using the first of Eq. (l23l) and exploiting the non-singularity of S2 
one has 

S2Va2eV^ S2 = S2lJ2nS2 =^ Va2nV'^ = (T2n, (137) 

which implies that V is a symplectic matrix (we are assuming a2i = o"2„). Moreover, 
from the second condition in Eqs. (1231 ) for S and S', we obtain 

S2as1W^ = S2Vas'[ , =^ s'^ = Ws4V, (138) 

because S2 is non-singular and V is symplectic. By considering the third condition 
(l23l) one then has 

Wissa2nsl + s,a2nsl)W^ = Wa2nW = a2n (139) 

which prove that PF is a symplectic. Finally, let us observe that the proof above does 
not use the non- singularity of S3. Indeed, one can relax this hypothesis and assume 
more simply that there exists a W such that S3 = Ws^; from Eqs. (1231 ) W has to still 
be a symplectic matrix but S3 and Sg may be singular. 

As an application of these equivalent unitary dilation results, we can find an al- 
ternative canonical form to the one in Sec. 12.51 with the same si and S4 but with S2 
and S3 of the following anti-diagonal block form 



Sj 



Fj 
G, 
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(140) 



where, for j = 2, 3, Fj, Gj are n x n real matrices. Imposing Eqs. (l23l) . one obtains 
the following relations 



J' 



F2G2 



F, - F^ 







J' 



IT 



FsGl 



1 



JG^ — G2J 



IT 



ni 

0, 



(141) 



the solution of which provides the following unitary dilation, 



S 



' In 





J 






G2 







.Gl 


-G^^tn 




-J) 


Un 






Gl^rG^^ 



(142) 



where again G2 is an arbitrary (non- singular) matrix and the eigenvalues of J are 
assumed to be different from 1. This solution is unitarily equivalent to the one in Eq. 
(|79l ) by applying V = —(y^n and 



W 



G^V-^Ga 

"Ct2 J Ct9 U 



(143) 



as above. 



E The ideal-like quantum channel 

Here we consider a quantum channel with X = t^n but F ^ with rank less than 
2n, which can be described in terms of only n additional (environmental) modes. We 
call it ideal-like quantum channel. Accordingly, the canonical unitary transformation 
IJ of Eq. (I20I) will be uniquely determined by a 4n x 4n real matrix S G Sp^An, R) 
of block form in Eq. (f22)) . where Sj are 2n x 2n real matrices. Particularly, si = S4 = 

l2n, 



S3 



F3 





G3 



S2 



-Gl 





-^3 



(144) 



with F3 and G3 being n x n real matrices such that F^^G^ = G^Fj, = 0, in order to 
satisfy the symplectic conditions in Eqs. (|23] ). Taking advantage of the freedom in 
the choice of the unitary dilation shown in Appendix |Dl the matrix S can be put in 



the form of Eq. (l22l) in which s[ = 54 = ! 



2n, 



4 





In 












(145) 



where F3 is assumed non-singular. In this respect, one uses V,W E Sp{2n, H) (of 
App. ini) of the following form 



V 



-F3 








(146) 



and W = V ^. Similarly, one can proceed, if G3 is non-singular, and obtains a 
similar structure for S as above. As concerns the weak-degradability properties, if 
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one assumes the initial environmental input state as 7^; = diag(2A^ + 1, 2M+ 1, 2A^+ 
1, 2M + 1), the eigenvalues of F - X^X-^{Y + ia)X-^X + ia are {2M, 2(M + 
1), 2N, 2(N + 1)}, which are always positive for any N ^ and M ^ 0; hence, this 
channel with 7^ as above is always weakly degradable. 

Finally, one may consider another ideal-like channel with X = 11 2n and Y = 
[(1 — (j3)/4]®", i.e. $x,y = <S>'i=ii^i)i^ where the single-mode Bi channel is defined 
in Ref. [ FTTl as X = I2 and F = (1 — o"3)/4. Trivially, this multi-mode channel is 
always WD (like Bi) and is able to transfer a quantum state without decoherence 
with the maximum quantum capacity (like for the single-mode case ifTTl '). 
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